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PHASE  I  FINAL  REPORT 
AESTRACT 


Concerning  adaptive  array  and  null  steering  applications  of 
polynomials  uiith  restricted  coefficients,  the  basic  mathematical 
question  to  consider  in  electronic  beam  steering,  with  a  discrete 
array  consisting  of  omn i-direct ional  elements  spaced  at  equal 
increments  along  a  straight  line,  is  how  coefficients  of  a  polynomial 
may  be  chosen  in  a  robust  yet  computationally  efficient  manner  so  as 
to  arrive  at  a  desired  beam  pattern.  In  numerous  applications,  these 
coefficients  are  required  to  satisfy  certain  restr ictions,  such  as  a 
bound  on  their  dynamic  range.  Thus,  particularly  in  null  steering,  it 
is  often  advantageous,  or  even  necessary,  for  the  shading  coefficients 
to  all  have  the  same  magnitude.  Basic  properties  of  such  polynomials 
and  their  applications  to  beamforming  are  described. 


1 . 


A.  STATEMENT  OF  WORK 


Determine  the  robustness  of  the  polynomials  with  unimodular 
coefficients  introduced  by  the  principal  investigator  in  1977. 

Determine  how  to  best  introduce  multiple  nulls  by  employing  the  theory 
of  polynomials  with  unimodular  coef f icients. 

Develop  adaptive  array  processing  algorithms  employing  the  concent  o- 
unimodular  shading  coefficients. 

Determine  to  what  extent  coherent  variation  of  shading  coefficients 
can  be  controlled  by  restricting  them  to  be  unimodular. 

Determine  methods  of  choosing  the  best  averaging  process,  as  described 
in  the  Phase  I  Work  Plan,  for  constructing  polynomials  with  the 
required  properties. 

Solve  the  Littlewood  and  Erdos  conjectures  for  F„. 

Construct  actual  polynomials  which  satisfy  the  Littlewood  conjecture 
for  G„. 

Develop  a  theory  of  interpolation  by  polynomials  with  unimodular 
coefficients  and  determine  applications  of  these  ideas  to  null 
steering  and  adaptive  array  processing. 

Study  the  zeroes  of  these  polynomials  for  both  theoretical  and  applied 
purposes. 


i  i  . 


3.  STATUS  OF  RESEARCH  EFFORT 


The  following  narrative  describes  the  significant  accomplish¬ 
ments,  progress  towards  research  objectives,  new  discoveries,  and 
specific  applications  of  the  Phase  I  research  carried  out  by 
Prometheus  Inc. 

I.  INTRODUCTION 


The 

fundamental  purpose  of  our 

current 

research  is 

to 

contribute 

to 

a  deeper 

understanding 

of 

the 

properties 

and 

appl ications 

of 

polynomials 

with  restricted 

coefficients,  and 

in 

particular 

with 

coefficients 

of  magnitude 

1. 

Such 

understanding 

1  s 

crucial  to  the  area  of  array  design  and  to  the  construction  of  robust, 
camputat ional ly  efficient  adaptive  array  algorithms.  This  report, 
with  its  appendices,  summarizes  our  progress  to  date. 

The  basic  mathematical  question  to  consider  in  electronic  beam 
steering,  with  a  discrete  array  consisting  of  omnidirectional  elements 
spaced  at  equal  increments  along  a  line,  is  how  coefficients  of  a 
polynomial  may  be  chosen  so  as  to  arrive  at  a  desired  beam  pattern. 
In  numerous  applications  these  coefficients  are  required  to  satisfy 
certain  restrictions,  such  as  a  bound  on  their  dynamic  range.  Here, 
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dynamic  range  refers  to  the  ratio  of  the  largest  to  the  smallest 
magnitude.  Thus,  particularly  in  null  steering,  it  is  often 
advantageous,  or  even  necessary,  for  the  shading  coefficients  to  all 
have  the  same  magnitude. 

Although  the  mathematical,  statistical,  and  physical  problems 
that  arise  in  the  consideration  of  array  shading  have  been  studied  for 
roughly  half  a  century,  many  interesting  questions  remain.  This  is 
true  even  for  the  ■simple1'  case  of  the  discrete  line  array  referred  to 
above.  In  this  case,  of  course  letting  n  denote  the  number  of 
elements,  the  pattern  function  G(z)  is  a  polynomial  of  degree  n-1,  2 
is  a  point  on  the  unit  circle,  and  the  shading  coefficients  are  just 
the  n  coefficients  of  this  polynomial.  An  important  reason  for 
performing  array  shading  is  to  shape  the  pattern  function  G  so  that  it 
has  low  sidelobes  and  small  beamwidth.  As  is  well  known,  both  of 
these  quantities  cannot  be  minimized  simultaneously,  and  the  choice  of 
shading  coefficients  results  in  a  tradeoff  between  these  two  desirable 
ends. 


Electronic  beam  steering  is  another  fundamental  purpose  of 
array  shading,  and  it  is  this  application  that  we  address  in  this 
report.  In  addition  to  permitting  the  rotation  of  the  main  response 
axis  of  the  pattern  function,  beam  steering  allows  the  simultaneous 
formation  of  a  number  of  beams  in  different  directions.  In 
particular,  if  sources  of  interference  lie  at  bearings  different  from 


that  of  the  desired  signal,  then  the  signal-to-noise  plus  interference 


ratio  (SNIR)  may  be  increased  dramatically  by  directing  nulls  of  the 
pattern  function  toward  these  interfering  sources,  in  spite  of  the 
fact  that  the  absolute  power  of  the  desired  signal  is  thereby  reduced. 
Adaptive  techniques  have  been  developed,  by  which  array  processing 
systems  can  electron i cal ly  respond  to  an  unknown  interference 
environment.  However,  although  the  basic  adaptive  array  principles 
have  been  known  for  some  time,  their  application  has  been  limited  by 
hardware  constraints  and  by  the  lack  of  sufficiently  robust,  real-time 
algorithms.  New  approaches  to  this  latter  consideration  are  described 
herein. 

There  are  many  cases  when  constraints  must  be  placed  upon  the 
magnitudes  of  the  coefficients  of  the  pattern  function.  Thus,  as 
explained  by  Hudson  C131,  when  coefficients  are  implemented  by 
attenuation  they  must  be  scaled  so  that  the  largest  modulus  is  unity, 
since  the  amplitude  gain  far  the  desired  transmission,  and  even  the 
overall  output  signal-to-noise  ratio  (SNR),  can  be  reduced  by  large 
coefficients.  In  discussing  main-lobe  constraints  on  optimal  arrays, 
Hudson  observes  that  when  a  main-lobe  null  is  created  very  large 
shading  coefficients  are  formed,  resulting  in  enhanced  output  of 
uncorrelated  noise.  Hence,  size  restrictions  on  the  coefficients  are 
again  required. 

On  the  other  hand,  in  a  situation  such  as  occurs  in  an 
adaptive  radar  receiver  after  clutter  has  decayed  due  to  increasing 
range,  so  that  there  will  be  few  and  widely  spaced  target  echoes  of 
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minimal  power  compared  to  a  steady  jamming  source,  it  is  necessary  to 


constrain  the  adaptive  array  so  that  the  shading  coefficients  are 
prevented  from  falling  to  zero.  A  similar  situation  occurs  in  an 
adaptive  antenna  using  the  least  mean  square  (LMS)  algorithm,  where 
the  shading  coefficients  will  decay  to  zero  if  either  the  signal  level 
falls  to  0,  or  if  the  reference  signal  is  absent  for  some  reason. 
One  method  of  controlling  this  is  to  substitute  the  steered  gradient 
system  described  by  Griffiths  far  the  reference  signal  LMS  antenna, 
but  this  has  the  disadvantage  of  being  very  sensitive  to  errors  in  the 


assumed  direction  of  the  desired  signal. 


As  mentioned  earlier,  another  approach  to  these  questions  is 
to  restrict  the  dynamic  range  of  the  shading  coefficients.  Although 
an  informal  rule  of  thumb  for  this  range  appears  to  be  *2  and  everyone 
is  happy,  10  and  some  are  happy,  100  and  nobody  is  happy, a  a  formal 
mathematical  study  of  the  relevant  properties  of  polynomials,  whose 
coefficients  are  thereby  restricted,  does  not  seem  to  have  been 
previously  undertaken.  An  important  thrust  of  the  research  effort 
reported  herein  has  been  to  initiate  such  a  study,  and  to  relate  the 
large  amount  of  work  that  has  been  accomplished  by  mathematicians,  on 
polynomials  with  restricted  coefficients,  to  the  above  applications. 
These  efforts  will  continue  m  Phase  II. 


Furthermore,  there  is  an  intimate  relationship  between  the 
engineering  questions  described  above  and  several  areas  of  classical 
mathematical  analysis.  Foremost  among  the  problems  of  mutual  interest 


is  tha  question  of  how  close  to  constant  the  modulus 


of  a  polynomial 


can  be  along  some  curve,  typically  the  unit  circle-  This  is  of  dees 


concern  to  theoretical  mathematical  analysts  because  of  the 


fundamental  nature  of  polynomials  and  the  simplicity  and  intrinsic 


beauty  of  the  question,  and  it  is  equally  important  to  engineers 


working  in  such  fields  as  -y  design,  adaptive  beamforming  and  null 


steering,  filter  design,  peak  power  limited  transmitting,  and  tha 


design  of  reflection  phase  gratings.  This  report  describes  cur 


research  into  both  aspects  of  this  remarkable  intertwining  between  the 


disciplines  of  pure  mathematics  and  engineering. 


I.  MATHEMATICAL  RESULTS 


Concerning  the  purely  mathematical  aspects  of  our  work,  note 
that  properties  of  polynomials  with  restricted  coefficients  have  been 
the  subject  of  much  fruitful  research  in  twentieth  century 
mathematical  aialysis.  Of  particular  interest  have  been  polynomials 
with  coefficients  +1  or  complex  of  modulus  one.  The  study  of  such 
functions  uias  apparently  initiated  by  G.H.  Hardy  (see  Zygmund  C23, 
p. 199  1 ) ,  and  furthered  by  J.E.  Littlewood,  P.  Erdos  and  others. 

For  the  purposes  of  this  discussion,  it  will  be  convenient  to 
introduce  the  notation  of  Littlewood  C17D.  Thus,  let  F„  and  G„  be, 
respectively,  the  class  of  ail  polynomials  of  the  form 

n  n 

f(z)=  5T  ±1  z  “ ,  g(z)  =  I  exp  (atai)  z*, 

t-o 

where fz^l  and  the  a*  are  arbitrary  real  constants.  Clearly  the  La  norm 
of  g  is  (n+l)^  far  all  g  4  Gn  (and  hence  for  all  f4  F„<G„),  and  the 
question  "how  close  can  such  a  g  come  to  satisfying 


The  first  qualitative  result  concerning  the  above  cuesticn  for 
Gn  was  obtained  by  G.H.  Hardy  C  23,  p.  199  3,  who  demonstrated  The 
existence  of  a  positive  constant  C  and  a  sequence  ^g„J  ,  ]n  ^  G„  , 
satisfying  ]g„(z)|  <  C  nw2  for  all  n  and  z.  The  identical  result  for 
F„  uas  obtained  by  Shapiro  C 22 3  and  published  by  Rudin  C19  3. 
Littleuood  C163  conjectured  that  there  exist  positive  constants  A  and 
E  such  that,  for  any  n,  there  is  an  f^F^Cg^G,,)  satisfying 

An  1  xa<  |f  (z)|  <Bn*'2  (An‘/a<  jg(z)|  <En  »  '2) 

for  all  z,  uihile  Erdos  conjectured  C  1 1  3  that  there  is  a  positive 
constant  C  such  that  for  n>2,  ||g|l  >  ( 1+C )  n 1  for  all  g<G„  (and  hence 

for  all  f<F„).  Analogous  conjectures  for  the  L*»  norms  of  g^5n  were 
settled  in  a  series  of  papers  by  Beller  and  Newman  (Cl  3,  C2  3,  and 
C33).  Beller  and  Neuman  C43  also  proved  the  Littleuood  conjecture  for 
polynomials  whose  coefficients  have  moduli  bounded  by  1,  after 
observing  that  the  proof  of  this  result  given  by  Clunie  CIO  3  depended 
upon  an  erroneous  result  of  Littleuood.  In  C153  Korner  was  able  to 
modify  the  results  of  Byrnes  in  C63  to  prove  the  Littleuood  conjecture 
for  G„,  and  then  Kahane  C143  showed  that  the  Erdos  conjecture  is  false 
for  G„.  These  conjectures  for  F„  remain  unresolved. 

One  approach  to  the  Erdos  conjecture  for  polynomials  in  F„  is 
to  consider  their  L*  norm.  Our  principal  results  in  this  regard  are 

described  in  Appendix  B.  For  polynomials  in  Gn,  we  have  the  following 

result : 


LbganeiD  !•  For  each  positive  integer  n,  there  is  a  sequence 


of  coefficients  £c„j  such  that  allfc*|  =  1  and 

k:0 


s# 

(  A 

^  r“ 


Z  c*  e1*^-  d*X<n  +  l  >a+4Cn+l  )='*. 

>uo 


PH2Qf  We  show  that  in  fact  the  Gauss  coefficients. 


satisfy  the  required  property.  Toward  that  end,  note  that 


r  n 

31  o„e 1  2!  cme  —  1 1 

fflso 


=  n  +  l  +  ZT  <  ST c«»j cm > e 1  J< 

j*<9  '•'to 


Therefore,  by  Parseval’s  Theorem,  assuming  for  convenience 


that  n  is  even. 


~9.tr  \  e  1  K,>|  dd-=(n  +  l)a+ 


jiO1  m;0 


1  C„ 


=  <  n+1 ) a  + 


a 


=<n+l )*  +  43 _ 


/  c.  uLZ 


i* '  \  JT/n  -ill 


<  ( n+ 1  )  a  +  4  ST  mini,*,  f  n  +  1  )  =*/ <  2  j  )  2 ) 
£<n  +  l )2  +  4 ( n+ 1 ) 3'2. 
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where  we  have  used  the  facts  that 


|sin  jx/sin  x|  <  j ,  and  ll/sin  xj  <  T7/2x  for  3<Tx<TT/2. 

This  completes  the  proof  of  Theorem  1. 

Another  method  of  constructing  polynomials  with  unimodular 
coefficients  is  to  form  a  suitable  weighted  average  of  existing  ones. 

For  example,  we  may  employ  a  slight  variation  of  the  basic 

construction  in  C61  as  follows: 

For  each  m,  0<m<Na-l,  and  for  z=eat_11*,  let 

#/-/.*>/  +  Q* j  +  fctJ 

P«(Z)=P„:0)  =  £  £  e  2  , 

Clearly,  each  P«.  is  a  polynomial  of  degree  Na-1  with  coefficients  of 
modulus  1.  Furthermore,  it  fallows  from  C63  that,  for  a  suitable  small 

positive  E  (i.e.,  of  order  N_a),  |P*,(0)|  is  essentially  flat  for 

E  -  m/Na  <  $  <  1 -£-m/Na  . 


Now  define  P*<€>)  by 

P-(6)  =  2_  z"'rJ 
m  zO 


P*  is  a  polynomial  of  degree  N*-l  with  coefficients  of  modulus  1, 
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Also,  by  writing 


P-(-e-) 


and  letting 


.  £ 


ZTT( 


;  A 


Ce 

*30  i-O 


"KT 


7J 


I  -  t 


-e. 


vtrl  A/  0 


=  -N-*< A+BN+CN*)  for  0<A<N~1,  0<B<N-1,  0<C<N2-1, 


it  is  seen  that 


P-<«  -  N«  " 


»  <-  rv 
n" 


so  that  |p*(-©)|  =  Na. 

In  addition,  the  essential  flatness  of  |P*(^>/  in  the  interval  ££<?-<l-E, 
where  now  EL  is  of  order  N-"*,  follows  as  before.  However,  numerical 
evidence  suggests  that  P~(  2-tN~*)  =0(1),  a  similar  situation  to  that 
which  occured  with  the  original  polynomials  [63.  This  being  the  case, 
P~  is  not  quite  a  Kahane-type  polynomial,  as  we  had  originally  hoped. 

Note,  however,  that  the  above  method  of  constructing  P*(-©-)  car. 
also  be  employed  to  create  new  flat  spectrum  sequences,  which  are 
periodic  sequences  E a„ 3^q  with  the  property  that  their  discrete 
Fourier  transform  (DFT)  has  a  power  spectrum  consisting  of  a  very 
small  number  (usually  1  or  2)  of  distinct  values.  This  is  because  the 
DFT  can  be  thought  of  as  the  values  of  the  polynomial 
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P(z)=  H  a  *2“, 
l:o 

where  n  is  the  period,  at  the  n-th  roots  of  unity.  Our  construct icn 
yields  polynomials  whose  spectra  are  essentially  flat  at  almost  all 
paints  of  the  unit  circle,  not  just  at  the  roots  of  unity.  Observe 
that  flat  spectrum  sequences  constructed  in  this  manner  satisfy  the 
additional  property  that  all  of  the  terms  of  the  original  sequence 
have  the  same  magnitude.  Applications  of  these  concepts  to  notch 
filtering  and  communications  are  discussed  elsewhere  in  this  report. 

Another  method  of  viewing  these  questions  is  in  the  context  of 
interpolation  problems.  As  noted  earlier,  for  any  P<G„,  the  Parseval 
Theorem  implies  that  the  La  norm  of  P  on  the  unit  circle  C  is 
(n+l)*'a.  Furthermore,  since  jP(z)|*  is  expressible  as  z-"Q(z),  where 
Q  is  of  deg  ee  2n,  there  can  be,  at  most,  2n  distinct  points  z*  where 

I P<  z* )  I  =  ( n+1 ) 


Let  us  call  such  a  set  of  points  an  La  Interpolating  Set  for  P.  A 
natural  question  is  which,  if  any,  subsets  of  C  consisting  of  2n 
points  can  be  an  L*  interpolating  set  for  some  P  of  the  required  form. 

In  its  full  generality,  this  question  appears  to  be  quite 
difficult.  For  n=l,  it  is  trivial  to  show  that  S=£a,b^  is  an  L2 
Interpolating  Set  if  and  only  if  b=-a.  For  arbitrary  n,  observe  that 


-11- 


TTTTT 


CTCTWWWTOTOgTTO A .V.WVWIW '■* m* -J- ^  ‘-'^^ -'^>^ 


3lV% 


for  S  =  $z*\  to  be  an  La  Interpolating  Set,  the  coefficients  of  P 

^  ->  K; o 


■Tiust  be  chosen  so  that 


!  (  s  )  —  <  n+1 )  zn  =  OC  fr  (=-z^)  , 


where  DC  is  a  constant  of  modulus  one.  Furthermore,  the  coefficient  of 
2"  on  the  left  side  of  this  equation  vanishes,  so  the  same  must  be 
true  on  the  right  side.  Clearly,  this  will  be  a  very  rare  occurence, 
so  that  most  sets  will  not  be  La  Interpolating  Sets.  In  fact,  it  is 
not  at  all  obvious  that  for  n>l,  there  exist  any  La  Interpolating 
Sets.  Thus  far,  we  are  only  able  to  show  that  if  S  is  to  be  such  a 
set,  its  elements  cannot  be  too  close  to  each  other.  More  precisely. 


Theorem  2  For  any  n,  there  is  an  E>0  such  that  no  S  of  the 


form 


S=Je^'}  ,  with  <  E.  for  l<k<2n. 


is  an  La  Interpolating  Set  for  any  P€G«. 


Proof  of  Theorem  2.  Assume  the  contrary.  Fix  n.  Then,  for 
any  E>0,  there  is  a  set 


S=S(E,n)  =  £e  ,  with|^»,|  <  E.  for  l<k<2n, 

such  that  S  is  an  La  Interpolating  Set  for  some  P,  say 
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Choose  a  sequence  of  positive  £.’s,  say^C-j^,  ,  approaching  0. 


For  each  k,  0<k<n,  the  sequence  Jag.,*) 

is  bounded,  and  all  terms  of  each  of  these  n+1  sequences  have  modulus 
1.  By  the  standard  method  of  choosing  a  convergent  subsequence  for  one 
k  at  a  time,  we  can  find  a  strictly  increasing  sequence  of  positive 
integers 


Jm_,  L,  ,  and  a  set  |aw^ 

of  complex  numbers  all  of  modulus  1,  such  that 

1 at  •  h \  converges  to  a*  for  every  k, 

t  ~i>  j-  , 

Since  )  ,  (eA^  )|  -(n+1 ) 1 -'a  can’t  change  sign  for 

|e^.|<exrr. 


0<k<n. 


uie  can  assume,  by  taking  another  subsequence  if  necessary,  that  either 


]  Pf  (e‘*  )  -  (n+ll,/J 


is  always  positive  or  always  negative  for 


e«v )  «?<rr. 


Suppose  the  former  (the  argument  being  the  same  in  the  latter  case), 
and  define 


P=(z>  =  Za„  z*. 


Clearly  (  z ) ^  converges  uniformly  to  Pa(z)  on  |z(  =  1, 

so  that  |Pc,(e‘®’)|  >(n+l  )  *'*  for  Q<&<2TT. 


Since  the  L3  norm  of  P«  is  (n+l)1'3,  this  is  impossible,  and  the  proof 
of  Theorem  2  is  complete. 


Also  of  interest  is  the  locations  of  the  zeroes  of  polynomials 

with  unimodular  coefficients.  This  is  directly  related  to  many  othen 

problems  discussed  herein  and  has  obvious  importance  in  the  choice  of 

pattern  functions  for  null  steering.  To  quantify  this  question,  let  r 
oC* 

ew  J,  l<j<n,  be  the  zeroes  of  P„  4-  G„,  normalize  P„  so  that  the 
coefficient  of  z"  is  1.  and  define 


where 


A 


max 


the  maximum 


min 

j 

is 


|  1  and 


max  ( 


r 

i-' 


k 

)  , 


ta^en  over  all  such  Pn(2). 
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Exnre  any  P  i  ( z  )=z-e  for  some  real  c<i.  .  it  is  obvious  that 


??=  =  0. 


Considering  the  case  n=2. 


lot,  I*.  L(**.**0 

Pa(z)=za— (r  i  +r2  e  >  z+r  tr  ae 


so  that 


r.ra=  rte  +rae 


***  j  =] 


Assume  that  rt>l.  Since 


*  1  ^ 

l=|r1e,‘  ’  +  rae11  1  >rl-r2  =  rL-  l/rt, 


the  maximum  value  for  r^l/r,  (hence  the  maximum  value  for  r 
achieved  when  r,-l/r!=l,  or 


r  1=  “Jt*  ,  and  r2= 


In  this  case, 


sTT-  \  J~s-\ 

r  t  - 1  =  — - —  and  l-r2=  * - ~=-  < 

1  v  fT 
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so  that 


x  U - 1  _  3- jr 
*’  jr-/  2 


Also, 


C 

/*"a.a  =  Max  (<r-l  )*-*-(  1-1/r  )a) 


By  an  elementary  calculus  argument,  it  is  seen  that  this  maximum 


occurs 


for  r  =  (  )  /2.  Thus, 


=  i  5-2,5  . 


We  leave  as  an  open  question  the  behavior  of  other  values  of  X,  and 

The  final  mathematical  question  which  we  examine  here  is  one 
which  is  particularly  important  for  applications.  Namely,  how  robust 
are  polynomials  with  unimodular  coefficients?  Thus,  expressing  P<G„  as 


P< z)  =  £  A *z«,  / Ate/ =  1, 

k.:0 


consider  that  the  coefficients  A*  are  not  fixed  quantities  but  may  be 
written  as  A*  -*■  a*«,  where  aw  =  a„(u i)  is  a  random  variable  representing 
the  coefficient  error.  Hence,  the  polynomial  becomes 


Q<z)=*  A„zw  +  Z.d 

k:0  kz° 
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inhere  the  final  sum  represents  an  error  polynomial  whose  behavior  is 
to  be  analyzed  to  determine  the  robustness  of  a  system  employing  these 
f unctions. 


Since  the  class  Fn  of  polynomials  of  the  form 


P(z>=  £  ±1: 
k'-O 


is  of  particular  interest,  let  us  begin  by  examining  error  polynomials 
of  the  form 


f  (  z  )  =  21  a^k 
kiO 


inhere  aM  are  real  random  variables  which  we  assume  to  be  independent 
and  normally  distributed  with  mean  0  and  variance  1,  to  be  written 
N ( 0, 1 ) .  Later,  we  shall  indicate  how  to  generalize  our  work  to  the 
case  of  other  random  coefficients  and,  finally,  to  complex  random 
coefficients. 


Let  us  rewrite  the  error  polynomial  as 


f<t)  =  a^lwle’ 
kzo 


where  the  ah  are  independent  real  random  variables  which  are  N(0,1) 
and  t  <  C0,2TI3.  We  begin  by  considering  the  real  part  of  f(t),  which 
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Re  f(t:=P(t,w)  =  £.a,«<w)cos  kt. 

k-o 


In  estimating  the  effects  of  coefficient  inaccuracy  in 
transversal  filtering,  Sersho,  Gopenath  and  Odlyzko  C121  analyze 
similar  error  polynomials.  In  comparing  various  performance  measures, 
they  conclude  that  the  most  appropriate  measure  is  the  maximum 
deviation  in  frequency  response  magnitude  from  the  desired  values  over 
a  particular  frequency  band.  Their  analyses  and  conclusions  apply  in 
our  case  and  lead  us  to  estimate  the  random  variable 


M*<w)  =  sup  P(t,w)  . 


Our  methods,  which  differ  from  those  of  Gersho  et  al,  are  based  on 
pioneering  work  in  random  polynomials  done  by  Salem  and  Zygmund  C201. 
They  yield  good  estimates  without  unknown  constants  and  are 
computationally  simple. 


Let  us  now  outline  the  procedure,  beginning  with  two  lemmas 
which  give  estimates  of  the  size  of  P<t,w). 


Lemma  JL.  E(e*•p,<*•'■',)  <  e 


-i-M 


,  where  E  represents  expectation. 


Proof.  It  is  well  known  that  for  a„~-N(0,l>,  the  moment 
generating  function  of  aw  is  given  by  E(e  )  =  e 


Then 


LtJk>k'A>J.<Ji*Ji>.li'jh<.l 


xP(t^)  c  D.akMc  oskt 

E <  e  J  )  =  EJe*-11 


r\  X  cos  kt  Q.(u>) 

=  TT  E(e  *  ) 


n  >a.  fu>)  cos  bt  , 
E(TT  e  J 

k=o 


X  (-os  l 

n  - g— 

TT  e  r 


n  i 
<  TTe  'L 

fc=0 


where  the  second  line  follows  by  the  independence  of  the  a^’s.  This 
completes  the  proof  of  Lemma  1. 


Lemma  2.  Define  M(w)  =  sup  |p(  t,w)/  . 

teloifl) 

Then  there  exists,  for  each  w,  a  random  set  I(w)  £  C0,2TT)  of  Lebesque 
measure  at  least  1/n  such  that 


P(t,w)|>  M(w)/2  for  each  te  I(w>. 


Proof.  By  Bernstein’s  theorem. 


Ip'II  =  sup  I  kf-  (t,w)l  <  nfl(w). 
00  ‘Ufyirfl  K  1 


Since  P(t,ui)is  continuous  for  each  w,  there  exists  t„cC0,2TT)  such 
that  M(w)  =  +P(t„,w).  Thus,  for  fixed  w. 


P(t,w)-P(ta,w)|  <  lip'll  <1  t-tofn  M(w). 


If  M(u>)  =  -*-P<t„,w>,  |P(t,w)-n(w)f  _P1(w)/2  for  /  t-t«Jn<l/2  (mod  2TT ) 


PS 


a  /.v.y.v.v.v.viv.v;  r.V.' »- 


If  M'w)  =  -»(t0,uj),  P  ( t ,  w )  +M  ( w )  <M(w>/2  for  |t-t„|n<l/2  (mod 


Thus,  Jp(  t,  w)|>M<w) /2  for  teC  tQ-l/2n,  t«»+l/2n]  mod  (2TT), 


and  Lemma  2  is  proven, 


We  combine  the  above  lemmas  to  obtain  the  principa 


estimation . 


\  (V/) 


Theorem  3.  E(eM%^2)  <4TTne 


Pr oo f .  E(eM^-'a)=  dP(u) 


J'  nl  i  (W) 


( t )  dt  dP(ui), 


where  we  have  used  Lemma  2  with 


li<«i(t)= 


i  if  t  (f  r (u j) 

o  ofUt  */tse. 


Therefore, 


T-ir 

E(e"^*)<nj  ^  j  e  «««-»?■'*  1  r  <«,  ( t  )dtj  dP(w) 


-n_  o 


r/  vr 

<n j  /  J  Csip<*'“'  +  e  «.*->  ]dtj  dP(  w) 


because 
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Now  use  Fubini’s  Theorem  "to  interchange  the  order  of 
integration,  and  apply  Lemma  1  to  obtain 


yr  %  c*rt) 

E(eM^a)<n  |  2e  dt=  4T7ne 

0 


and  we  have  Theorem  3. 


To  obtain  our  estimate  of  the  size  of  M,  we  prove  th  = 
:o i lowing : 


Theorem  4. 


-«x/?  (n+i) 

P  £M>  °<  ]  <  4TIne  for  o<  >  0, 


Let  ^->0  and  use  the  Markov  Inequality  and  Theorem  3  to 


obtain 


P  EM  >  =  P(e  "X'*  >  e*%/S  <  E(e"*'*> /e**^  <  4TTne 


%-M-  T 


To  get  the  best  bound  for  a  given  o< ,  we  employ  elementary  calculus  to 
minimize 


h  <P-)  =  Xa<n+1  )  /2-  a«</2, 

yielding  7-e=  °< /2 (n+1  ) .  Substituting  ^  =  into  the  right  side  above 
yields  Thecrem  4. 
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We  can  non)  obtain  similar  estimates  for 


N ( u )  =  sup  |f(t,w)|  ,  where 

teiozn) 


f<t ,  id)  =  £  ak(u))e  1(<* 


and  the  ak  are  independent  N(0,1). 


Theorem  5. 


P^N(w)  ><<|< 


STTne 


r  n 

fn SQl-  Let  Re  f  =  £  a*  cos  kt  and  Im  f  =  £a„  sin  kt.  We  h< 

k'-°  k* 

shown  in  Theorem  4  that 


Pj  sup  /Re  f(t,w)|  >  <*  /J2f  <  4TTn 

tt)  i 


-  «1/ 1 6  (n-h  i) 


Ey  applying  the  analysis  in  Lemmas  1  and  2  and  Theorems  3  and  4  to  Ir 


f,  uie  obtain 


PVsup  / Im  f(t,w)/  >  ®<  /J~2  1  <  4TTi 


-°r7/£  frW) 


Since  |f  ( t ,  w )  j  2  =  / Re  f(t,w)ja  +  / 1 m  f(t,w)|a,  we  have 

£?i >  °<a3c(supC|Re  f(t,w)|  2  >«C2/2jUEsup  /im  f(t,w)| 


Theref ore, 


>  «<2/23. 


STTne 


p(Nk  •</  <  Fjf  sup  |Re  f(t,u)|  >  «  /|T2/+  p/ sup  llm  f  Ct,u)/  >  «:  /  lz] 

V/,*Ou,) 


.  <  _  r.*,  •>  -  v  •  ,  *'•  *  _  •>  .  «  _» 


§ 


completing  the  proof  of  Theorem  5. 

We  now  choose  K  to  obtain  suitable  explicit  bounds.  If  we  put 
o<=  4  J  C  n  +  1 (  log  STTn+4 loc  log  n  )  .  then 


pf N  >  <  1  /  ( log  n  )  •*, 


which  improves  an  estimate  in  Theorem  1  of  £121. 
If  o<  =  4i(n+l)  log  (8TTna),  then 


’[N  >  <J< 


In  this  way,  estimates  of  any  desired  precision  may  be  obtained. 

When  coefficient  error  is  due  to  digital  round-off,  it  is 

reasonable  to  assume  that  the  coefficients  a*  are  uniformly 

distributed  on  the  interval  C -1/2,  1/21,  denoted  ^ -1  /2,  1  /2  1.  The 

fallowing  result  gives  the  analogue  of  Theorem  5  in  this  case. 

Theorem  6.  Let 
r\ 

f(t,w)  =  2La„(w)e  “•*,  t  e  £0,217), 

k-o 

where  the  a*  are  independent  random  variables  which  are  UC-1/2,1/23. 


Then 


><j< 


8TTne 


Proof.  A  direct  calculation  shows  that 


E(e  )  = 


e.  -  -  e- 


r  )/%  if  ^  ^ 


1  *  5  =o 


By  simple  use  of  the  exponential  series,  we  have 


using  2**k  !  <  (2k)'  (2k+l) 


Finally, 


>ak 

E  ( a  ) 


We  now  fellow  the  steps  in  the  procedure  for  norma 
coefficients  to  obtain  similar  results  to  Theorem  5.  First  consider 


P(t,w>=  2L  a*  cos  kt  and  M(w)  =  sup  |P(t,w)|  . 

ku>  tejoiv) 

Fallowing  Lemma  1,  we  obtain 

(*+l) 

E(ex-* *•“* )  <  e 


and  from  Lemma  2, 

z. 

2jr  (n+i) 

E(e  )  <  4TTne 


As  in  Theorem  4,  we  have 


r  <x/l(mi) 

P]M  >  «<;<  4TTne 


We  obtain  a  similar  bound  for  the  maximum  a*  the  imaginary  part  o 
f(t,w),  and  we  combine  these  as  in  Theorem  5  to  obtain 


f  /  >  <K 


SM 

8TIne 


! 

i 


I 

i 
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which  completes  the  proof  of  Theorem  6. 

Remark.  We  may  obtain  specific  estimates  for  P^N  >  in  the 
uniform  case  just  as  we  did  for  the  normal  case  after  Theorem  5.  For 
example,  setting 


o<  =*  2  J<  n  +  1 ) log  ( 3TTna) ,  we  obtain  P(N  >*0  <  1/n. 

We  are  now  prepared  to  give  estimates  on  the  error  polynomials 
obtained  from  polynomials  with  unimodular  coefficients.  These  error 
polynomials  are  of  the  form 

n 

f(t,w)  =  Z  an  e  **•*, 

k*> 

where  the  ak  are  independent  complex  random  variables.  In  other 

words,  a«.=<«-*-ij#,.  where  o<*  and  are  independent  real  random 

variables.  We  state  the  main  Theorem. 

Theorem  7.  Let  N(w)=sup  /f(t,w)/,  where 

f(t,w)=  X  awt(w)e1*‘* 
to 

and  the  a „  =»<h.+i>^w  are  independent  random  variables.  We  suppose  that 
the  real  and  imaginary  components  of  each  a*,  namely  «<«  and  are 

independent.  Then 


a)  If  <Ku  and  are  N<0,1),  P^N  >  *<]  <  16TTne 


(b)  If  «<„  and  A,  are  U  C -1/2,  1/21,  p{ N  >  o<  ]  <  l^TTm 
Proof .  Let  us  write 

f(t,w)  =  f i(t,w)  -*■  ifa(t,w>  ,  where 


-o fV/6  (nil) 


f  i  (t,w>  =  £.<><*  eiw* 
fc=o 


and  fa(t,w>  =  Z.  A.e1**. 
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Consider  the  case  ithere  A<  are  independent  N(0,1). 

N t  ( u» )  =  sup  ’  If  i  ( t,  ui)|  ,  i  — 1,2,  then 
ttfrjtn) 

(if  ( t,u»>|  >  c  (if  t  (t,ui)/  >  °<  /  2>u//fa(t,ui>/  >  o<  /  2/  , 
since  I  f /  <  |f  »  |  +  |f  a  / .  Therefore, 

p£.N(w)  >«K^<  P^Nt(w>  >  «<  /  2^+  P^Naiw)  >  *  /  2^  =  SP^Ntdu)  >  o<  /  2^ . 


Applying  Theorem  6  yields 


P^N(ui)  >  <X j  <  16TIn 


uhich  is  (a).  A  similar  procedure  applied  to  the  uniform  case  gives 
(b),  and  uie  have  Theorem  7. 

In  summary,  we  have  obtained  easily  calculated  and  flexible 
upper  bounds  on  the  maximum  deviation  caused  by  random  error 
polynomials.  These  bounds  offer  estimates  of  the  robustness  of  a 
discrete  array  of  omnidirectional  elements. 


III.  APPLICATIONS 


As  mentioned  in  the  introduction,  applications  of  polynomials 
with  restricted  coefficients  abound  in  the  engineering  jjor  Id.  Those 
wnich  we  focus  on  herein  include  null  steering,  adaptive  beamfcrmmg, 
notch  filtering,  peak  power  limited  transmitting,  and  the  synthesis  cf 
low  peak-factor  signals  and  flat  spectrum  sequences. 

In  Appendix  C,  several  new  designs  of  analytic  null  steering 
algorithms  for  linear  arrays  are  described.  Two  of  them,  the  £— 
Technique  and  the  Positive  Coefficient  Model,  allow  for  placing  an 
arbitrary  number  cf  nulls  in  arbitrary  directions,  while  maintaining 
main  beam  and  sidelobe  level  control.  A  method  of  incorporating  these 
deterministic  null  steering  techniques  into  existing  adaptive 
algorithms  is  proposed.  The  resulting  Direct  Adaptive  Nulling  System 
offers  the  possibility  of  significant  increases  in  array  performance 
at  very  little  cost.  This  possibility  will  be  investigated  in  full 
detail  in  Phase  II. 

A  major  reason  for  combining  deterministic  methods  with 
existing  techniques  is  that  arrays  must  ordinarily  deal  with 
significant  random  noise.  In  these  cases,  one  has  no  a  priori 
information  about  the  direction  or  nature  of  such  unwanted  signals. 
Thus,  in  such  applications,  as  well  as  in  cases  where  advance 
knowledge  of  jammer  character i st i cs  is  lacking,  indirect  statistical 
methods  are  unavoidable,  although  their  efficiency  may  be  greatly 
increased  by  combining  them  with  analytic  approaches. 

There  exist  applications,  however,  where  much  is  known  in 
advance  about  the  character  1 st 1 cs  of  both  the  desired  signals  and  the 


undesired  noise.  This  is  especially  true  where  one  has  control  of  the 
generation  of  these  waveforms.  Thus,  in  the  case  where  one  system  is 
producing  both  offensive  signals  (i.e.,  searching  for  and  homing  in  on 
targets)  and  defensive  signals  (i.e.,  identifying  and  tracking 
incoming  weapons),  sc  that  mutual  interference  becomes  a  predominant 
concern,  the  problem  is  almost  exclusively  determin i st i c  in  nature.  In 
such  cases,  robust  and  computationally  efficient  analytic  algorithms 
controlling  both  the  individual  performance  of  the  offensive  and 
defensive  signals  and  the  interactive  jamming  between  them  are  crucial 
to  mission  success.  The  application  to  such  cases  of  the  new 
determin ist i c  null  steering  algorithms  developed  by  Prometheus  will  be 
analyzed  in  detail  during  Phase  II. 

A  related  problem  is  the  determination  of  optimal  shading 
coefficients  for  a  conformal  array.  As  is  well  known,  using  various 
measures  of  optimality,  this  is  a  computational  problem  of  order  n3, 
where  n  is  the  number  of  array  elements.  Thus,  the  computational  load 
will  be  reduced  by  a  factor  of  8  if  the  coefficients  may  be  restricted 
to  be  real.  Circumstances  where  this  occurs  are  described  in  Appendix 
D.  A  different  method  of  improving  computational  efficiency,  namely  a 
convex  programming  approach,  will  be  an  important  focus  of  our  Phase 
I  I  research. 

Another  interesting  application  of  our  concepts  is  to  notch 
filters.  Appendix  E  describes  a  nearly  ideal  notch  filter  employing 
coefficients  of  equal  magnitude.  Applications  to  the  design  of 
transmitting  antenna  arrays  are  discussed  briefly.  The  construction 
is  based  upon  earlier  worn,  of  the  author  involving  polynomials  with 
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restricted  coefficients.  The  fundamental  idea  employed  in  Appendix  E 
to  construct  a  notch  filter  with  a  single  notch  may  be  combined  with 
the  concept  of  an  n-nomial  C5  3  tc  produce  nearly  ideal  filters  with 
multiple  notches.  Furthermore,  as  noted  elsewhere,  sera  coefficients 
do  not  affect  the  dynamic  range,  so  that  these  multi-notch  fi Iters 
maintain  the  property  of  having  unit  dynamic  range.  Details  of  this 
new  construction  and  experimental  results  employing  it  will  be 
developed  during  Phase  II. 

In  addition  to  their  use  in  the  construction  of  notch  filters, 
Byrnes  Polynomials  16,14,153  have  potential  applications  to  the  design 
of  peak  power  limited  transmitters  and  the  synthesis  of  low  peak- 
factor  signals  and  flat  spectrum  sequences.  In  transmitter  design,  for 
example,  one  is  often  faced  with  a  peak  power  constraint.  Under 
various  conditions,  the  transmitter  output  may  be  modeled  as  a 
polynomial.  Here  the  maximum  modulus  of  the  polynomial  on  the  unit 
circle  represents  the  peak  power,  while  the  L2  norm  of  the  polynomial 
is  the  average  power.  Thus  the  classical  engineering  problem  of 
minimizing  the  peak  to  average  ratio  becomes  the  mathematical  question 
of  minimizing  the  ratio  of  the  sup  norm  to  the  L2  norm  of  a  polynomial 
on  the  unit  circle. 

In  the  trivial  case  where  one  frequency  is  to  be  transmitted 
(ie,  the  polynomial  can  be  a  mononomial),  clearly  the  ideal  value  1 
for  the  peak  to  average  ratio  is  achieved,  and  the  polynomial  is 
indeed  of  constant  modulus  on  the  unit  circle.  For  the  more 
interesting  and  practical  case  of  transmitting  many  linearly 
increasing  frequencies,  it  is  usually  desired  to  transmit  each 
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frequency  at  the  same  power,  which  should  be  as  large  as  passible.  As 
the  power  of  each  individual  frequency  is  represented  by  the  modulus 
of  the  corresponding  coefficient,  the  mathematical  question  naturally 
arises  of  how  close  to  constant  the  modulus  of  a  polynomial  with 
equimodular  coefficients  can  be  on  the  unit  circle. 

More  precisely,  if  n  pure  tones  are  transmitted  with 
frequencies  of  the  form  f0+k4,  where  f0  is  the  fundamental  frequency 
and  A  is  the  increment,  then  the  waveform  is 
x(t)  =  21  A„  cos  C  2FT  (  f  o+H)  t  +dk  3 

fco 

=  |S(t)|  cosCarg  S(t)  +  2TT  f „t  1. 

Here,  S(t)  =  Z.  A*e  etarl,*A*,^i«  =  phase  and  Ate=  power  in  k**  tone. 

<-o 

As  mentioned,  almost  always  all  frequencies  are  transmitted 
with  equal  power,  so  that  Ata=l.  To  minimize  the  peak  power  of  x(t), 
the  maximum  (over  t)  of  /x  <  t  >  /  must  be  minimized  (over  O,. ) .  It  is 
relatively  straightforward  to  see  that  the  exact  problem  is  to  obtain 

I  n-l  tfl  _  / 

min  maxi  JF  e  e‘arfk^t  , 

^  t  1  t*  ' 

a  job  which  is  performed  by  the  Byrnes  polynomials  C51  in  nearly  ideal 
fashion . 

The  adaptation  of  such  polynomials  to  these  problems  is 
important,  since  in  applications  like  the  Link  11  Communications 
System,  the  average  power  is  usually  maintained  at  one  tenth  or  less 
of  its  theoretical  ideal  to  prevent  transmitter  overload.  Employing 
concepts  such  as  those  described  above  should  yield  a  significant 
reduction  in  the  peak-to-aver age  ratio,  thereby  allowing  a  large 
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increase  in  average  power,  hence  a  more  efficient  communications 
system.  These  considerations  also  show  that  the  Byrnes  construct  ion 
has  direct  application  to  the  synthesis  of  low  peak-factor  signals. 

Now  consider  the  problem  of  designing  a  flat  spectrum  sequence 

JO 

{ak}  ,  as  defined  on  oage  13.  These  sequences  have  direct  use  in  such 
diverse  areas  as  concert  hall  acoustics,  the  quieting  of  an  object’s 
response  to  radar  and  active  sonar,  and  speech  synthesis.  Schroeder 
C211  presents  many  of  the  fascinating  details  of  these  applications. 

As  we  observed  earlier,  the  DFT  can  be  thought  of  as  the 
values  of  the  polynomial 

n- 1 

P(z)=  X  a.:*, 

where  n  is  the  period,  at  the  n-th  roots  of  unity.  The  Eyrnes 
construction  C  6  1  yields  polynomials  whose  spectra  are  essentially  flat 
at  almost  all  points  of  the  unit  circle,  not  just  at  the  roots  of 
unity.  Furthermore,  they  have  the  additional  property  that  all  of  the 
terms  of  the  original  sequence,  ^a*^  have  the  same  magnitude. 
Applications  of  these  concepts  to  notch  filtering  and  communications 
are  discussed  elsewhere  in  this  report. 

In  our  final  application,  we  have  begun  to  exploit  the  great 
success  of  J.P.  Kahane  C141  in  salving  the  Littlewood  conjecture.  As 
we  note  in  part  II,  Kahane  showed  that  there  indeed  exist  polynomials 
with  unimodular  coefficients  whose  modulus  is  essentially  constant  on 
the  unit  circle.  It  is  our  opinion  that  the  breakthrough  of  Kahane 
was  due  to  the  ingenious  use  of  randomness  and  probability  in  his 


construct  ion .  Behind  his  and  previous  approaches  was  the  idea  of 
Gauss,  viz.  the  "Gauss  Sums."  To  put  it  quite  simply,  we  feel  that 
Littlewood’s  problem  mas  vanquished  by  the  "equation" 

Kahane  =  Gauss  Sums  +  Probabilistic  Choices. 

Our  idea  is  to  exploit  the  Kahane  breakthrough  by  developing 
methods  to  judiciously  make  the  " Probabi 1 istic  Choices"  referred  to 
above,  and  thereby  convert  Kahane’ s  "randomized"  proof  into  a 
constructive  one.  This  mould  not  only  result  in  exciting  new 
mathematics,  but  mould  also  be  directly  applicable  to  several 
important  engineering  problems.  In  addition  to  the  areas  of  peak 
power  limited  transmitting  and  flat  spectrum  sequences  discussed 
earlier,  such  polynomials  mould  find  immediate  use  in  the  design  of 
reflection  phase  gratings,  and  therefore  be  employable  in  solving 
concert  hall  acoustics  problems  and  in  quieting  the  response  of  an 
object  to  sonar  or  radar.  Another  potential  application  of  this 
"educated  randomness"  construction  is  in  the  synthesis  of  multielement 
omnidirectional  beam  patterns. 

In  the  concert  hall  acoustics  application  of  reflection  phase 
gratings,  it  is  desired  to  design  the  ceiling  so  that  sound  is  widely 
scattered  except  in  the  specular  direction.  As  described  earlier  and 
in  Appendix  E,  in  the  context  of  notch  filter  design,  the  Byrnes 
polynomials  C63  place  a  null  in  any  given  direction  while  the 
coefficients  maintain  their  other  desirable  properties  of  being  both 
flat  spectrum  and  lorn  correlation  sequences.  Thus,  they  might  even  be 
preferable  to  the  Kahane  polynomials  in  this  context.  This  also 
appears  to  apply  to  monostatic  radar,  where  the  null  mould  be  placed 


in  the  direction  of  the  radar.  For  bistatic  radar,  on  the  other  hand, 
the  receiver  direction  is  often  unknown.  Thus,  if  a  construction 
based  upon  the  Kahane  polynomials  could  be  employed,  radar  energy 
would  be  reflected  equally  in  all  directions,  theredy  reducing  the 
probability  that  there  would  be  enough  energy  reflected  in  any 
particular  direction  to  enable  detection.  A  possible  undersea 
application  of  these  ideas  occurs  in  the  design  of  baffles  used  to 
quiet  machinery  noise  from  submarines,  in  an  attempt  to  prevent  the 
noise  from  escaping  the  hull.  Note  that  our  constructions  would 
complement  the  coatings  that  are  already  in  use,  or  being  designed,  to 
attack  these  problems,  since  these  coatings  provide  uniform 
attenuation.  Furthermore  surface  structures  based  upon  the  Byrnes 
polynomials  would  have  the  highly  diffusing  property  over  a  large  set 
of  frequencies.  It  is  not  yet  clear  whether  the  Kahane  polynomials 
also  yield  this  important  property.  The  design  of  two-dimensional 
arrays  so  that  energy  may  be  scattered  with  equal  intensity  over  the 
solid  angle  is  also  of  considerable  interest.  It  appears  that  a 
straightforward  product  formulation  gives  the  desired  results  for  the 
Byrnes  polynomials,  but  the  situation  is  not  so  clear  for  Kahane 
polynomials.  Our  Pnase  II  research  will  focus  upon  the  many 
fascinating  questions  raised  in  this  final  paragraph. 
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APPENDIX  A 

Additional  Information 
A- I.  TECHNICAL  PUBLICATIONS 

1.  Eyrnes,  J.S.,  The  Minima:-;  Optimization  of  an  Antenna  Array 
Employing  Restricted  Coefficients,  Scientia  (to  appear,  1937). 

2.  Newman,  D.J.  and  Byrnes,  J.S.,  The  L*  Norm  of  a  Polynomial  with 
Coefficients  +1,  (submitted). 

3.  Byrnes,  J.3.  and  Newman,  D.J.,  Null  Steering  Employing  Polynomials 
with  Restricted  Coef f icients,  (submitted). 

4.  Byrnes,  J.S.,  A  Notched  Filter  Employing  Coefficients  of  Equal 
Magnitude,  (submitted). 

A— I I .  PROFESSIONAL  PERSONNEL 

1.  James  S.  Byrnes,  Principal  Investigator.  President,  Prometheus  Inc. 
Ph.D.  (Mathematics),  Yeshiva  University,  1967. 

2.  Donald  J.  Newman,  Principal  Scientist.  Ph.D.  (Mathematics),  Harvard 
University,  1953. 

3.  Stephen  Boyd,  Senior  Scientist.  Ph.D.  (Electrical  Engineering), 
University  of  California  at  Berkeley,  1985. 

4.  Andre  Giroux,  Senior  Scientist.  Ph.D.  (Mathematics),  University  of 
Montreal,  1973. 

5.  Martin  Goldstein,  Senior  Scientist.  Ph.D.  (Mathematics),  University 
of  Wisconsin,  1969. 
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A- III.  INTERACTIONS 

I.  Invited  Papers  Presented  by  J.  3.  Byrnes. 

1.  Department  of  Electrical  Engineering,  Stanford  University. 
7  October,  1936. 

2.  Naval  Research  Laboratory,  Washington,  DC,  12  Novemcer 
i9e6. 


3.  Naval  Underwater  Systems  Center,  New  London,  CT. .  28 

January  1987. 

II.  Upcoming  Invited  Lectures  on  this  Work,  by  J.  S.  Byrnes. 

1.  Department  of  Mathematics,  University  of  Maryland,  April, 
1987. 

2.  NATO  Advanced  Study  Institute  on  EM  Modeling  ard 

Measurements  for  Analysis  and  Synthesis,  Tuscany,  Italy. 

August,  1987. 

3.  XXIInd  General  Assembly  of  the  International  Union  of  Radio 
Science,  Tel  Aviv,  Israel,  September,  1987. 

III.  Consultations  on  Potential  Air  Force  and  Navy  Applications  of 

Prometheus  Ideas. 

1.  Lecture  presented  by  Donald  J.  Newman  at  Raytheon  Equipment 

Division,  Uayland,  MA,  26  September  1986.  Also,  Newman  and 

Byrnes  of  Prometheus  consulted  with  Duane  Matthiesen,  Hero 
Groginsky,  Leon  Wardle,  and  Fred  Daum  of  Raytheon,  Hans 

Steyskal  and  John  K.  (Jay)  Schindler  of  RADC  -  Hanscomb,  Len 
Smith  of  MITRE,  and  Richard  Turyn  of  GTE  -  Natick. 

2.  Stanford  University,  7  October  1986.  Byrnes  and  Bey d 

consulted  with  Professor  Bernard  Widrow. 

3.  Naval  Research  Laboratory,  10  November  1986.  Byrnes 
consulted  with  Emanuel  Vegh,  Paul  Crepeau,  Frank  Kretschmer , 
Karl  Gerlach,  J.  Rao,  Jack  Ahern,  S.N.  Samaddar,  and  Allen 
Mi  1 ler . 

4.  Various  Locations,  15-17  December  1986.  See  attached  Trip 
Report . 

5.  Naval  Underwater  Systems  Center,  28  January  1987.  Byrnes 
consulted  with  Roger  Dwyer,  John  Fay,  Kurt  Hafner,  Huck  Qua:; . 
Ed  Eby,  and  Cliff  Carter. 
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Prometheus  Inc. 

103  MANSFIELD  STREET,  SHARON,  MA  02067 


401-849-3389 
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617-784-2355  1-800-225-5040 


TRIP  REPORT 

In  connaction  with  contract  #  F4962Z-86-C-0088 
NULL  STEERING  APPLICATIONS  OF  POLYNOMIALS 
WITH  UNIMODULAR  COEFFICIENTS 
Jamas  S.  Byrnes,  Principal  Investigator 


A  Phasa  I  SBIR  Contract  with  the 
Air  Force  Office  of  Scientific  Research 
Or.  Arje  Nachman,  Program  Manager 


Oates  of  Trip:  13-17  December,  1986 

Prometheus  Personnel:  James  S.  Byrnes,  Principal  Investigator 

Donald  J.  Newman,  Principal  Scientist 
Stephen  Boyd,  Senior  Scientist 

Installations  Visited: 

1.  Rome  Air  Development  Center,  Hanscomb  Air  Force  Base,  MA  (Byrnes, 
Newman,  Boyd). 

2.  MIT  Lincoln  Labs,  Bedford,  MA  (Byrnes,  Newman,  Boyd). 

3.  Raytheon  Equipment  Division,  Uayland,  MA  (Byrnes,  Newman). 

4.  MITRE,  Bedford,  MA  (Byrnes,  Newman). 

3.  Raytheon  Submarine  Signal  Division,  Portsmouth,  RI  (Byrnes). 

Purposes  of  trip:  To  discuss  new  results  obtained  by  Prometheus  under 
the  above  contract,  to  learn  of  passible  applications  of  these  results 
at  the  above  installations,  and  to  obtain  suggestions  for  directions 
of  future  research. 


Personnel  Visited: 

1.  RADC:  Dr.  Robert  Mailloux  (primary  point  of  contact).  Dr.  Hans 
Steyskal  (617-377-2032),  Dr.  Robert  Shore,  Mr.  Jeff  Herd. 

2.  MIT  Lincoln  Labs:  Dr.  Charlas  Rader  (617-863-5300,  *2374). 

3.  Raytheon-Way land:  Dr.  Eli  Brookner  (primary  point  of  contact, 
617-338-2721,  x3636),  Dr.  James  Mullen,  Mr.  Fred  Daum. 

4.  MITRE:  Dr.  Dean  Carhoun  (primary  point  of  contact,  617-271-2318), 
Dr.  Irving  Reed,  Dr.  Warren  Wilson  (617-271—3913),  Dr.  John  Cozzens 
(617-271-3404),  Mr.  Len  Smith  (617-271-3905),  several  others. 

3.  Raytheon-Portsmouth:  Dr.  Dave  DeFanti  (primary  point  of  contact, 
401-847-8000,  x4411 ) ,  Dr.  Stan  Chamberlain,  Dr.  Roger  Pridham,  Mr.  A1 
Serheim. 

Results:  Most  of  the  above  individuals  found  the  Prometheus  work  very 
interesting,  and  had  many  suggestions  for  applications  and  follow-on 
work.  These  suggestions  have  been  incorporated  into  our  proposal 
"Polynomials 'with  Restricted  Coefficients  and  their  Applications,"  and 
will  form  an  integral  part  of  our  Phase  II  proposal. 
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APPENDIX  3 

The  L'*  Norm  of  a  Polynomial  With  Coeff icients  ±1  ' 

Donald  J.  Newman  and  J.  S.  Byrnes 
Prometheus  Inc. 

103  Mansfield  Street 
Sharon,  MA  02067 

A  classic  unresolved  question  regarding  n— th  degree 
polynomials  with  coefficients  ±1  is  whether  the  maximum  modulus  of 
such  a  polynomial  on  the  unit  circle  can  be  n  ‘■/*+o( n  *'a) .  As  shown 
by  Kahane  C23,  if  complex  coefficients  of  modulus  1  are  allowed  then 
not  only  is  it  possible  for  this  property  to  be  satisfied,  but  the 
minimum  modulus  can  be  nwa+o(nwa)  as  well.  Specifically,  Kahane 
proved  that  for  any  n  there  is  a  polynomial  of  degree  n  with 
coefficients  of  modulus  one  whose  modulus  everywhere  on  the  unit 
circle  is  n*'atrO(na,''**log  n). 

Erdos  C 1 1  had  conjectured  the  existence  of  a  c>0  such  that, 
for  any  polynomial  P  of  the  types  described,  I|p|/  >(l+c)n‘/a.  Clearly 
the  Kahane  result  disproved  this  conjecture  for  the  modulus  1  case, 
but  the  situation  for  coefficients  ±1  remains  open.  Employing  an 
elegant  construction  Shapiro  C4,33  demonstrated  the  ach ievabi 1 i ty  of 
the  order  of  magnitude  nwa,  but  the  maximum  modulus  of  the  Shapiro 
polynomials  is  (2n)*'a.  Motivated  by  these  considerations  we  examine 
the  L*  norm  of  such  polynomials.  As  one  might  expect,  this  leads  to 
several  interesting  combinatorial  questions.  We  provide  answers  to 
some  of  these,  and  conclude  with  a  refined  version  of  the  Erdos 


conjecture 


rw.N^TV'.v-V'.v^rjvvv'.’T'^rjv', 


Throughout  the  paper  n  will  be  a  positive  integer,  P(z)  will 
denote  a  polynomial  of  degree  n-1  with  coefficients  +1,  and  z  will 
lie  on  the  unit  circle.  Thus, 


r  k,  ,  zm& 

P(z)=^£^z  ,  each  =1  or  -1,  z=e  ,  0<€K1. 


All  integrals  will  be  over  €><C0,11.  We  begin  with  a  Lemma. 


T  r - 

Lemma  II  Pjl  -  )  £  i  £  t 

n  V  j  k  L  " 

^  IIPIIM'  iwnfj* 

/hr/  /t/hri 


;  nr  V‘)a\<k)(iv-%^n 

y  j ^  J  6  'Wo  ' 

*  ^constant  terms. 

Since  a  constant  term  occurs  in  this  product  if  and  only  if  j+k=l+m, 
the  result  fallows  immediately. 

4 

Of  interest  is  the  expected  value  E(  |/p||  >,  if  the 

£ 

coefficients  £j  are  chosen  at  random. 

4  f 

Theorem  1  E(  l|Pj|^)=2na— n. 

Proof  Clearly  if  exactly  3  of  the  indices  j,k,l,m  are 
identical,  or  if  at  least  3  of  them  are  different,  then  E(£ 

f  ,  J 

1 )=0.  It  therefore  follows  from  the  lemma  that 
L  "> 

Of  k-t 

For  each  of  the  n(n-l)/2  pairs  of  integers  p,q,  0<p<q<n, 


there  are  4  terms  appearing  in  (1),  namely  ^ ^  *7*  ’ 

,  Th*  only  other  terms  in  (1)  are  £*,  0<p<n.  Since  all 

of  these  terms  equal  1, 


E (  II pH  )=n+4(n )  ( n-1 ) /2=2n2-n, 

L* 

completing  the  proof  of  Theorem  1. 

•  We  now  observe  the  improvement  that  is  achieved  when  this 

random  choice  is  replaced  by  the  Shapiro  coefficients. 

Theorem  2  If  n=2“  and  P(z>  is  the  Shapiro  polynomial  of 
degree  n-1,  then 


||p||  »(4na-<-l )  “n  ) /3. 

a 


Proof  Shapiro*3  polynomials  are  defined,  together  with  his 
auxiliary  polynomials  Q,  by  the  recurrence  formulas 
P«<  z )=Q.<  z )  =  1 , 


P^x(z)=P»(  z  )+za  Qn.(  z  ) , 

(  z )  =*P«(  z  ) Qm(  z  ) ,  m>0. 


As  a  result. 


so  that,  as  is  well  known. 


*  _mtl 
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Now  (2>  and  (3)  yield 


'-KJl 


Next  we  observe 


that  z  Cl  P  is  composed  solely  of 


frequencies  which  are  positive  powers  of  z,  so  that  it  can  be  thought 
of  as  Q  P  %  where  P  is  the  "reversed*  polynomial  of  P  .  Thus 

Since  fsQP  is  analytic  and  0  at  the  origin, 

o* 

so  that  1 

f(fkf)z  --  t  $ISl'---tfiQj'lFjx-tflPj'  IqJ. 


Altogether  then,  we  have 


fun'* 


l"+3  -iflPJ*. 


The  remainder  of  the  proof  is  now  simply  induction  on  k.  The 
result  is  obviously  true  for  k*0,  since  Po<z)»l.  Furthermore,  from 


(5)  and  the  inductive  hypothesis 

SlPJ  --  ~ — j1  ~'z 


it  follows  that 


ft/- 


xk+1  z 
2  -  Z  - 


^7  (kw)t  Z 
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This  completes  the  proof  of  Theorem  2. 


Note  that  Theorem  2  implies  that  the  L*  norm  of  the  n-lm* 
degree  Shapiro  polynomial  is  asymptotic  to  v/rT  times  the  fourth  root 
of  4/3*1  •  0743  7  i/tT.  Based  upon  extensive  numerical  evidence  employing 
the  Bose*~Einstein  statistics  methodology  of  Statistical  Mechanics,  we 
conjecture  that  the  Shapiro  polynomials  do  not  give  the  minimum  L- 
norm  among  all  polynomials  of  the  same  degree  with  coefficients  ±l, 
but  that  this  minimum  L*  norm  is  asymptotical  ly  v'rT  times  the  fourtn 


root  of  6/5r>»l . 04664^'n'.  Observe  that  the  truth  of  this  conjecture 


would  imply  that  of  the  Erdos  conjecture  mentioned  earlier,  with 
c— ( 6/5  )  1 04664. 
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ABSTRACT 

Several  neu  designs  of  analytic  null  steering  algorithms  for 
linear  arrays  are  described.  Tuo  of  them,  the  /-Technique  and  the 
Positive  Coefficient  Model,  allow  for  placing  an  arbitrary  number  of 
nulls  in  arbitrary  directions,  while  maintaining  main  beam  and 
sidelobe  level  control.  A  method  of  incorporating  these  deterministic 
null  steering  techniques  into  existing  adaptive  algorithms  is 
proposed.  The  resulting  Direct  Adaptive  Nulling  System  offers  the 
possibility  of  significant  increases  in  array  performance  at  very 
little  cost. 


-  *■«  »  .  /.  »  .  r.  J*  «/* 


WK***!ry*fm*!*!  h .  l  .  .  mm  .  .  .  mi  s  » ■  j  j  j 


■w y'-.'i'.1  '..’V’T,*1..11  V  WVWWfrlYV*  "J1  ’J*  V  V-V.,,y.’yry7WA'  irrr:i,:T’.'»VFVwv»i 


Null  Steering  Employing  Polynomials  with  Restricted  Coefficients 
J.  S.  Byrnes  (Sr.  Member,  IEEE)  and  Donald  J.  Newman 

Prometheus  Inc. 

I.  INTRODUCTION 

In  view  of  the  well-known  one-to-one  cor respcndence  between 
polynomials  and  linear  arrays  with  commensurable  separations  between 
elements,  as  described  in  detail  by  Schelkunoff  C  9  ],  null  steering 
questions  involving  such  arrays  translate  directly  into  mathematical 
problems  regarding  the  locations  of  zeroes,  on  the  unit  circle,  of 
polynomials.  Furthermore,  physical  and  electronic  limitations  placed 
upon  the  array  elements,  such  as  a  maximum  allowable  power  or  a  bound 
on  the  dynamic  range,  imply  various  restrictions  upon  the  coefficients 
of  these  polynomials.  Here,  dynamic  range  refers  to  the  ratio  of  the 
magnitudes  of  the  largest  to  the  smallest  weight,  or  shading 
coefficient,  of  the  array.  Thus,  the  theoretically  challenging 
question  of  the  placement  of  zeroes  at  specified  points  on  the  unit 
circle,  of  polynomials  whose  coefficients  satisfy  certain 
restrictions,  is  also  a  problem  of  strong  practical  interest  to 
antenna  designers. 

The  design  of  filters  is  another  application  in  which  such 
questions  arise.  For  example,  the  classical  mathematical  problem  in 
natch  filter  design  is  to  produce  a  polynomial  whose  magnitude  on  the 
unit  circle  is  close  to  constant  in  almost  all  directions,  but  which 
has  a  small  number  (:e,  1,  2  or  3)  of  deep  nulls  ("notches")  at 
specified  points.  In  15]  the  construction  of  C4 3  is  employed  to 
produce  such  a  polynomial  having  one  null,  with  the  added  feature  that 
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all  coefficients  have  the  same  magnitude.  Hence,  the  dynamic  range  of 
the  notch  filters  presented  in  C53  is  one. 

This  paper  addresses  the  null  steering  application  describee 
in  the  first  paragraph.  There  are  several  factors  which  must  be 
considered  in  the  design  of  null  steering  algorithms.  In  addition  tc 
the  basic  problem  of  placing  the  nulls  the  main  beam  must  be  steered, 
the  width  of  the  main  lobe  controlled,  and  the  sidelobe  levels  must  be 
sufficiently  below  that  of  the  main  lobe.  Control  of  the  sidelobe 
level  is  usually  achieved  by  attenuating  the  shading  coefficients  as 
one  moves  away  from  the  center  of  the  array.  Often  these  attenuation 

factors  (Chebyshev,  Taylor,  etc. >  are  chosen  in  advance,  and  may  not 

be  easily  altered  once  the  array  is  in  place.  This  leads  directly  to 
a  beautiful  mathematical  question,  similar  to  the  peak-factor  problem 
in  engineering  attacked  by  Boyd  C2  3,  Schroeder  C103  and  others: 

Given  the  magnitude  of  the  coefficients  of  a  polynomial  P,  a 
finite  subset  S  of  the  unit  circle  C,  and  a  point  p<C  distinct 
from  those  in  S,  choose  the  phases  of  these  coefficients  so 

that  P<z)=0  for  all  z<S,  the  maximum  on  C  of  iP(z)i  occurs  at 

z=p,  and  the  maximum  of  ]p(z)i  on  a  subset  of  C  excluding  an 
appropriate  interval  (the  beamwidth)  around  p  is  as  small  as 
possible. 

We  consider  various  suboroblems  in  this  paper.  Research  on  the 
general  question  is  continuing. 

II.  DIRECT  ADAPTIVE  NULLING 

Currently  the  most  widely  used  class  of  null  steering  methods 
1  s  k  noun  as  adapt  1  ve  nulling  C'.  ,3,6,7,8,121.  Adapti  ve  arrays  have 
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developed  over  the  past  twenty-five  years  as  the  oref erred  method  of 
reducing  the  performance  deter icrat ion  in  signal  reception  systems 
which  is  inevitably  caused  by  uncesired  noise  entering  the  system. 
Sources  for  this  noise  include  multipath  affects,  electror.ic 
countermeasures,  clutter  scatterer  returns,  antenna  location  errors, 
array  element  thermal  noise,  etc.  The  prol i f eration  of  such  ncise 
sources  has  greatly  increased  the  importance  of  interference 
suppression  in  essentially  all  applications.  Although  such  adaptive 
methods  as  the  Uidrow  least  mean  squares  (LMS)  and  Howel ls-Applebaum 
sidelobe  canceller  have  achieved  considerable  success,  difficult 
problems  remain.  Foremost  among  these  are  poor  transient  response, 
signal  cancellation  resulting  from  interaction  between  signal  and 
interf erence,  excessive  computation  time,  and  sidelobe  degradation 
when  jammer  cancellation  is  attempted.  A  secondary  problem  is  the 
lack  of  control  in  adaptive  algorithms  of  the  dynamic  range  of  the 
weights. 

These  methods  are  indirect  adaptive  schemes^  they  do  nor 
explicitly  form  an  estimate  of  the  directions  of  arrival  of 
interfering  sources  or  explicitly  steer  nulls  in  those  directions.  A 
scheme  in  which  these  two  tasks  are  actually  performed  can  be  called  a 
direct  adaptive  algorithm.  Thus,  one  approach  to  the  solution  of  suer, 
problems  is  to  complement  an  appropriate  indirect  adaptive  algorithm 
with  the  analytic  null  steering  methods  described  herein.  In  this  way, 
the  actual  noise  suppression  achieved  can  be  enhanced  beyond  that 
which  would  be  available  through  either  adaptive  or  analytic  methods 
exclusively. 
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The  first  step  in  this  “Direct  Adaptive  Nulling  System"  wouia 


be  to  employ  available  techniques,  such  as  the  maximum  entropy  method 


spectral  estimation,  a  saatial  discrete  Fourier  transform,  of  the  arrai 


outputs,  a  search  in  angle  with  an  auxiliary  beam  C31,  or  aha 


loop  of  the  indirect  adaptive  algorithm,  to  estimate  noise  and  signal 


directions.  Our  analytic  methods  would  then  be  applied  to  choose 


shading  coefficients  which  place  nulls  in  the  estimated  noise 


directions,  while  maintaining  other  desired  properties  of  the  tea' 


pattern.  Feeding  these  coefficients  back  into  the  adaptive  portion  01 


the  algorithm  then  results  in  reestimates  of  the  directions  of  th- 


noises  and  signal,  which  are  used  in  turn  by  the  analytic  portion  to 


.mprove  the  choice  of  the  shading  coefficients.  This  process 


continues  until  convergence  is  achieved.  Furthermore,  each  execution 


of  the  analytic  step  is  essentially  instantaneous,  as  the  evaluation 


of  the  coefficients  given  the  required  directions  is  simply  a  matter 


of  plugging  the  data  into  elementary  formulas. 


We  expect  such  a  direct  scheme  to  outperform  the  indirect 


methods  in  cases  where  noise  or  interference  is  highly  correlated  witt 


the  desired  signal,  such  as  occurs  with  multipaths  or  intelligent 


Jamming.  Indirect  schemes  tend  to  perform  poorly  in 


environments.  In  addition, a  direct  scheme  allows  much  greater  use  o; 


prior  knowledge,  such  as  known  jammer  locations  or  known  multipaths. 


Thus,  this  interactive  direct  method  offers  the  possibility  if 


significant  increases  in  performance  at  very  little  cost. 


III.  COEFFICIENTS  OF  EQUAL  MAGNITUDE 


An  important  subproblem  of  the  general  mathematical  quest  icr, 
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described  earlier  is  the  case  when  all  of  the  coefficients  of  the 
polynomial  have  the  same  magnitude,  which,  by  normalization,  we  can 
assume  to  oe  one.  Such  phase-only  shading  occurs,  for  examp_e.  m  the 
design  of  transmitting  arrays  which  are  omn i di r act i cna 1  except  for 
specified  nulls.  These  features  are  crucial  in  certain  communications 


areas,  where  it  is  desired  to  null  out  listeners  in  known  directions 
while,  at  the  same  time,  for  maximum  efficiency,  all  antenna  elements 
are  broadcasting  at  full  power.  Also,  in  order  to  minimize  the 
relative  sice  of  the  quantisation  steps  in  a  gradient  algorithm  such 
as  Lh!S,  the  coefficient  magnitudes  should  be  kept  as  close  as  possible 
to  unity  C7,  p.  1531.  Note  that  this  * equimagn i tude“  property  of  the 
coefficients  precludes  the  use  of  attenuators,  with  the  concomitant 
savings  in  electronic  hardware. 

The  most  elementary  example  of  the  above  is  the  unshaded  array 
-  all  coefficients  are  1.  In  spite  of  its  simplicity,  this  uniform 
array  is  of  practical  importance.  Observe  that  in  this  case  the 
zeroes  of  the  polynomial  are  almost  uniformly  spaced  around  the  unit 
circle,  occurring  at  all  of  the  n+1-*  (where  n  is  the  degree  of  the 
polynomial)  roots  of  unity  except  z=l ,  where  there  is  a  maximum. 

At  the  other  extreme  is  the  case  wnere  an  n-fold  zero  is 
required  at  one  point.  One  application  of  this,  as  discussed  by 
Steyskal  till,  is  to  broaden  a  pattern  null  so  as  to  null  an  entire 
sector.  Clearly,  by  a  simple  change  of  variables,  this  zero-point  can 
be  assumed  to  occur  at  z=l.  It  is  a  straightforward  matter  to 
construct  such  a  polynomial  with  coefficients  of  magnitude  1;  in  fact, 
the  coefficients  may  all  be  taken  to  be  +1.  Namely,  define  P(z)  by: 
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P  (  Z  )*  =  T  I  (  1~Z3  )  (1) 

m=0 

The  problem  with  this  construction  is  that,  although  P(  z )  obviously 
satisfies  the  required  properties,  it  does  so  at  very  high  cost. 
Since  P  has  degree  2°-l,  its  realization  requires  an  array  with  2" 
elements.  We  show  in  the  following  theorem  that,  for  all  but  small 
values  of  n,  this  situation  may  be  greatly  improved  by  allowing 
coef  f  i  cier,  ts  to  be  0  as  well  as  +1.  Since  this  simply  means  chat  some 
array  elements  are  turned  off,  the  dynamic  ranee  of  the  coefficients 
is  not  affected  in  any  meaningful  wau. 

Theorem  X  Let  n>10.  Then  there  is  a  polynomial  Pic),  c“ 
degree  less  than  n3,  such  that  P(z)  has  an  n-fold  zero  at  z= 1 ,  and  all 
the  coefficients  of  P  are  either  +1  or  0. 

Proof  of  Xiisarem  X  Given  n>10,  choose  k  so  that 
2“>f<"2  (2). 

Since  x/ln(:<)  is  increasing  for  x>e,  and  10/ In  (  10  )  >3/ In  ( 2 )  ,  the  choice 
k  =n 3  certainly  implies  (2).  Actually,  for  any  03,  it  is  clear  that 
for  n  large  enough  we  may  take  k=na”“*,  but  such  precision  is  net 
necessary  here.  Now,  for  each  arbitrary  subset  S  of  the  set  of 
nonnegativs  integers  less  than  k,  let 

0  <  z  >  =  £L  cm, 

*|  eS 
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and  form  the  vector 


(Q<1>,  Q  ’  (  1 ) ,  Q 7  7  (  1  )  /  2  :  _ _ Q«"-‘  ’  (  1  )/  (n-1  )  1  . 


These  are  integer  vectors,  and  the  largest  entry  is  bounded  Dy 


k-l 

mn_ 


1  <  k  " . 


Thus, there  are  less  than  kn  such  vectors.  Since  there  are  2“  subsets 
S,  and  so  2“  polynomials  Q,  equation  (2)  implies  that  at  least  too 
distinct  polynomials,  say  Qt(z)  and  Qa(z),  have  the  same  associated 
vector.  Hence, P ( z ) =Q i ( z ) -Ga ( z )  is  the  desired  polynomial,  and  Theorem 
I  is  proven . 

The  idea  underlying  equation  (1),  which  we  call 


•encapsulation,"  may  also  be  employed  to  construct  polynomials  witr, 
coefficients  of  magnitude  1  that  place  any  number  of  arbitrary  nulls. 
Namely,  we  have: 

n 

Theorem  IX  For  any  positive  integer  n,  let  £zm3  be  an 

*>:/ 

arbitrary  set  of  (not  necessarily  distinct)  points  on  the  unit  circle. 
Then  there  is  a  polynomial  P(z>  with  coefficients  all  of  magnitude  1, 


of  degree  2"-l,  satisfying  P<zm)=0,  l<m<n. 

Pr oof  of  Theorem  11  As  indicated  above,  we  simply  produce  an 
explicit  formula  for  P(r): 


n  -l  ■» 

p<  z  >=  rr  i 2 


w- 1 

■Z-*  ) 


It  is  straightf oruard  to  see  that  this  P(z)  satisfies  the  require; 
properties,  establish ing  Theorem  II. 

IV.  SIMULTANEOUS  NULL  STEERING  AND  MAIN  BEAM  PLACEMENT 

In  the  previous  section^  uie  attacked  the  suborcblem  of  the 
general  question  stated  earlier  which  arises  when  the  dynamic  ranee 
of  the  coefficients  is  required  to  be  one.  We  now  consider  another 
aspect  of  the  original  problem,  perhaps  of  more  interest  to  antenna 
designers.  Namely,  how  can  arbitrary  nulls  be  placed  while 
maintaining  a  specified  main  beam  direction  and  specified  maximum 
sidelobe  level?  We  describe  two  methods,  the  /^-Technique  and  twe 
Positive  Coefficient  Model,  of  achieving  these  goals. 

To  set  the  problem,  again  let  n  denote  any  positive  integer, 
n 

and  let  S=Cz*3  be  an  arbitrary  set  of  (not  necessarily  distinct; 
paints  on  the  unit  circle.  Also,  let  zm  be  a  point  on  the  unit  circle 
distinct  from  those  in  S.  Our  methods  allow  the  placement  of  zeroes 
of  a  ploynomial  P  at  all  paints  in  S,  while  simultaneously  having  the 
maximum  of  1 P i  on  the  unit  circle  occur  at  z=z».  Furthermore,  the 
difference  between  |p<z«)i  and  the  highest  sidelobe  can  be  mace 
arbitrarily  large.  As  will  be  seen,  the  costs  encountered  i  ■- 
achieving  the  last  property  are  an  increase  in  the  degree  of  P,  anc  a 
loss  of  control  of  the  dynamic  range  of  the  coefficients. 

To  proceed  with  the  constructions,  define  the  angles  C«hm3. 
0_:m<n,  by  z»=exp (  idm )  ,  -TI<0*,<1T.  As  before,  a  simple  change  c* 

vanables  allows  us  to  assume  0«=0,  so  that  z«=l. 

Method  The  ^-Technique 

n 

Let  A=-cat~ 1 1  2  cot  (^m/2  )  3,  z*=exp(  ii), 

n-.  i 


C-9 


and  defire  Q(z)  by 


n 

Q ( z )  =  ( s-:*)TT  : . 

"1  s/ 

A  straightf oruiard  calculation  shows  that  Iq(z)1  has  a  relative  maximum 
at  z  =  l.  Hence,  for  c  a  large  enough  positive  integer", 
P  (  z  )  =  (  1  +  z  )  *ST*Q  (  2  )  will  certainly  satisfy  the  required  properties.  It 
can  be  shown  that,  in  order  to  guarantee  an  absolute  maximum  of  lp(z)| 
at  z=l,  it  is  sufficient  to  take  where  <=min iomi.  Of  course, 

in  order  to  further  increase  the  main  lobe  level  relative  tc  the 
sidelobes,  it  will  be  necessary  to  take  c  larger. 

Method  ,11.  The  Positive  Coefficient  Model 

Far  each  m,  l<m<n,  choose  the  smallest  positive  integer  k* 
such  that  exp(ikm^m)  lies  in  the  left  half  plane,  and  define  P(z)  by 


Clearly  P(z)  has  the  necessary  zeroes.  Furthermore,  all  of  the 
coefficients  of  P  are  positive,  so  that  the  maximum  of  |p<z)1  on  t-.e 
unit  circle  obviously  occurs  at  z=l.  Once  again, it  is  simple  to 
further  increase  the  main  lobe  level  relative  to  the  sidelobes  by 
multiplying  PCz)  by  an  appropriate  positive  integer  power  of  1+z. 

There  are  two  additional  points  which  can  be  made  about  the 
Positive  Coefficient  Model.  One  is  that  its  electronic  implementation 
will  be  greatly  simplified  as  compared  to  that  of  arbitrary  shading 
coefficients,  since  the  positivity  eliminates  the  need  for  phase 
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xrVXVT irnv^ji'vv* jv  vv  Al’vm ■vw  . v  v  v 


shifters.  A  second  is  that  some  control  of  the  dynamic  range  of  the 
coefficients  can  be  achieved  by  combining  this  method  mitt,  the 
encapsulation  technique  discussed  earlier,  if  me  again  ignore  the 


effects  of  3  coef f i cien ts . 


V.  CONCLUSION 


Various  aspects  of  a  fascinating  problem  in  classical 
mathematical  analysis,  with  direct  applications  to  antenna  array 
design,  have  been  discussed,  and  several  results  obtained.  Foremost 
among  these  are  two  analytic  methods  for  placing  an  arbitrary  number 
of  nulls  in  arbitrary  directions,  while  maintaining  .main  bean  and 
sidelobe  level  control.  A  method  of  incorporating  these  analytic  null 
steering  techniques  into  existing  adaptive  algorithms  is  proposed. 
The  resulting  Direct  Adaptive  Nulling  System  offers  the  possibility  of 
significant  increases  in  array  performance  at  very  little  cost. 


W.V  V  V 
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It  Is  well  known  that  the  determination  of  optimal  shading  coefficients 
for  an  antenna  array,  using  various  measures  of  optimality,  Is  a  computational 
problem  of  order  n^,  where  n  is  the  number  of  array  elements.  Consequently, 

If  It  Is  possible  to  achieve  optimum  results  by  employing  real  shading 
coefficients,  so  that  the  effective  number  of  coefficients  to  calculate  is 
halved,  the  computational  load  can  be  reduced  by  a  factor  of  8.  It  Is  shown 
by  Lewis  and  Strelt  [1]  that  real  coefficients  do  Indeed  suffice  for  the 
minimax  design  of  a  linear  antenna  array.  The  purpose  of  this  note  is  to 
prove  that  a  similar  result  can  be  achieved,  under  certain  circumstances,  for 
a  conformal  array. 

To  precisely  define  the  problem,  assume  that  there  are  n  elements  located 
at  points  (xj,  yj,  zj).  l<J<n,  and  a  set  F  of  farfield  points 
u-(cosa,  cosD,  cosy!  cos2a  ♦  cos20  +  cos2y  -  1.  Let  the  shading 
vector  be  w«(wi ,W2, . . .wn) ,  and  denote  the  beam  pattern  by 


T(u.w) 


n 

£ 

j-i 


exp 


(XjCOSo 


y^cosB 


ZjCOSY) 


(1) 


where  x  Is  the  wavelength  of  the  design  frequency.  The  minimax  optimization 
problem  can  be  stated  as  follows: 


(*)  Choose  shading  vector  w 
subject  to  the  normalization 


so  as  to  minimize  the  quantity 


Max 


T(u,w) 


T(ug.w)-1. 


(2) 


Here,  ug-(cosag,  cosOn,  cosYg)  is  a  fixed  farfield  point.  The 
normalization  (2)  Is  intended  to  f.orce  the  Maximum  Response  Axis  (MRA)  to 
occur  at  the  farfield  point  ug.  Consequently  ug  should  not  lie  In  the  set  F. 


i 


i 

i 

* 


It  Is  proven  herein  that  in  certain  cases  problem  (*)  can  be  solved,  with 
the  same  minimum  achieved,  when  the  weight  vector  w  Is  restricted  to  be  real. 
Namely,  the  fbllowlng  theorem  holds: 


•A 


'.UsY . 


Theorem  J.  Let  F  be  symmetric  with  respect  to  the  origin  of  (x,y,z)  space; 
that  Is,  u«f  If  and  only  if  -u*F.  Let  the  positive  z-axls  be  the  MRA,  and 
assume  that  all  Zj  are  Integer  multiples  of  x/2.  Then  the  minimum  achieved 
In  solving  problem  (*),  when  the  weight  vector  w  Is  restricted  to  be  real.  Is 
equal  to  the  minimum  achieved  without  this  restriction. 

Proof.  Let  v-(vi,V2 . vn)  be  any  (complex.  In  general)  vector 

which  solves  problem  (*),  and  let  vj«r<  ♦  Isj.  The  proof  will  be 

finished  If  It  can  be  shown  that  the  (real)  vector  r»(ri,r2 . rn)  also 

solves  (*) .  Toward  this  end,  first  consider  the  normalization  constraint 
(2).  Since  the  MRA  Is  the  positive  z-axls. 


oq  »  8q  ■  */2  and  yq  ■  0* 


Combining  this  with  (1)  and  (2)  yields 


1-T(u0,v)-^  Vj  exp  (1*J  Zj), 

J-l 


(3) 


so  that 


Z  ri  cos(?  'j)-  !j  *in(K)-  '• 


(4) 


Z  rj  s,n(^  sj  cos(^  zj)-  »• 


However,  since  all  zj  are  a  multiple  of  x/2.  It  follows  that  all  terms  In 
(4)  Involving  s1n(2wZj/x)  are  zero.  When  combined  with  (3)  this  Immediately 
Implies  that 


T<vr) -E  rj exp  ( ^  zj)  -  ^ 

j-i 


(5) 


which  means  that  the  vector  r  satisfies  the  normalization  constraint  also. 
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0  • 


To  see  that  r  satisfies  the  required  minimax  property,  recall  that  F  Is 
symmetric  with  respect  to  the  origin,  and  apply  the  method  of  [1].  Thus, 
letting  an  overbar  denote  complex  conjugate,  any  weight  vector  w  must  satisfy 


Max  |T(u,w) 

-  Max  |T(-u,w)| 

-  Max  T(u,w)j 

UcF  | 

ucF  |  1 

ucF  | 

so  that 

Max  JT(u,Re(w))|  *  Max  |T(u,  Iw  +  rw)! 
ucF  I  I  ueF  I  6  c  I 

<  \  Max  |t(u,w)|  +  ^  Max  |t(u,w)|  ■  Max  |T(u,w) 
ucF  *  ucF  ucF  ' 


Applying  (6)  to  the  optimum  vector  v  yields 

Max  T(u,r)j  <  Max  |t(u,v)| 
ucF  1  *  UcF  '  I 


However,  since  v  minimizes  this  maximum  among  all  normalized  weight  vectors, 
and  because  r  satisfies  (2),  It  follows  that  the  Inequality  In  (7)  must  In 
fact  be  an  equation.  Thus  the  real  vector  r  solves  problem  (*),  and  the 
theorem  Is  proven. 
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Me  now  show  that  the  at) ove  method  also  yields  a  more  general 
result  for  a  linear  array  than  that  in  C 1 3.  Namely,  we  prove  that  the 
shading  coefficients  in  a  minimax  optimisation  problem  for  a  linear 
array  can  be  taken  to  lie  along  any  fixed  line  L  through  the  origin  in 
the  complex  plane,  if  the  normalization  condition  is  altered  to  require 
that  the  sum  of  the  weights  be  either  of  the  two  fixed  complex  numbers 
of  modulus  one  lying  on  L.  Thus,  consider  n  omnidirectional  elements 
located  at  arbitrary  fixed  positions  xfc  along  the  x— axis.  The  minima;: 
problem  for  a  linear  array  steered  broadside  becomes: 

(*■*)  Choose  the  complex  weight  vector  w3(w(  ,  wx  ,  . .  .  ,  w  >  so  as  to 
minimize  the  quantity 


As  indicated  above,  Lewis  and  Streit  show  that  problem  (*«-)  has 
a  solution  with  a  real  weight  vector  w.  Our  general ization  of  their 
result  is  as  follows: 


Theorem  2.  Replace  the  normalization  (8)  by 

r> 

T(0)=  £  w  =exp(  ii*) .  (9) 

k;l  * 

Then  problem  (•**)  has  a  solution  where  all  wk  are  taken  along  the  line 
given  (in  polar  coordinates)  by 

Proof.  The  projection  of  each  component  wfc  of  any  weight  vector 
w  upon  the  line  is  given  by 

exp<  i'f )  ( cos  +  Re  w  +  sin  'f  Im  w  > , 

P  * 

so  that  the  method  employed  to  obtain  (7)  yields 
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<J-Max|^w  exp(-idfc  u|  + -i- Max|  £  w^  exp(-id^  u)[  =Ma>:|  £  w  e:<p(-xd^u| 

^  U  t-i'*'  U  kri  ^  k-l  K 


Furthermore,  it  follows  from  (9)  that 


r  i 

Lexp  ( i't')  ( cost  Re  w  +sin  f  In  u  )=exp(i't')  as  well. 

t  k 

krl 

Now  suppose  that  [v^]  is  any  sequence  of  weighting  coefficients 
satisfying  this  new  constrained  minimax  problem.  That  is, 


and,  for  any 


^  v  =exp <  ) 

r  k'-‘  k 

Vui  S  satisfying 

iui  f  . 

L  uk =exp( iT> , 


n  k:i 
Max  I  I  v  e 

K  v .  k 


I  I  '  I 

„  v  exp(-id  u)  <  Max|  w  exp(-id  u)| 
1/  k  k  ~  1  k-l  k  k 


Then,  by  the  above, 


exp<  i 'i')  <  cos  f  Re  v  +  sin  f  In  v  ) 

k  k 


also  satisfies  the  problem,  and  these  coefficients  all  lie  on  the  line 
d=T-  This  completes  the  proof  of  Theorem  2. 

Finally,  as  Lewis  and  Streit  note,  under  some  circumstances  a 
condition  such  as  T<u)=0  for  some  point  u*0  might  be  a  crucial 
requirement  added  to  problem  (**).  They  correctly  observe  that  now  a 
solution  with  real  coefficients  need  not  necessarily  exist.  However,  it 
is  quite  possible  that  a  judicious  choice  of  'f  in  the  above  will  allow 
for  the  solution  of  this  modified  problem  to  be  found.  We  leave  this  as 
a  subject  for  future  research. 
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Abstract 


A  nearly  ideal  notch  filter,  employing  coefficients  of  equal 
magnitude,  is  described.  Applications  to  the  design  of  transmitting 


antenna  arrays  are  discussed  briefly.  The  construction  is  based 


upon 


earlier  work  of  the  author  involving  polynomials  with  restricted 
coefficients. 


I.  INTRODUCTION 


The  classical  mathematical  problem  in  notch  filter  design  is 
to  produce  a  polynomial  whose  magnitude  on  the  unit  circle  is  close  to 
constant  in  almost  all  directions,  but  which  has  a  small  number  (i.e., 
1,  2  or  3)  of  deep  nulls  ('notches’)  at  specified  points.  Such 
filters  are  applied,  for  example,  to  remove  spectral  lines  from 
otherwise  broadband  spectra.  In  this  paper,  we  produce  such 
polynomials  having  one  null,  with  the  added  feature  that  all 


coefficients  have  th« 


same  magnitude.  For  convenience,  this  magnitude 


is  assumed  to  be  one.  Observe  that  this  ’unimodular’  property  allows 


a 


*„  *.  *.  *.  V-  ’v  n 


the  direct  application  of  these  polynomials  to  the  design  of 
transmitting  antenna  arrays  which  are  omnidirectional  except  for  1 
null.  This  feature  is  crucial  in  certain  communications  areas,  where 
it  is  desired  to  null  out  one  listener  in  a  known  direction  while  at 
the  same  time,  for  maximum  efficiency,  all  antenna  elements  are 
broadcasting  at  full  power. 


If  the  polynomial  P(z)  is  of  degree  n-1,  it  is  clear  from  the 
Parseval  Theorem  that  its  L2  norm  (i.e.,  RMS  value)  is  exactly  nw2, 
since  there  are  n  coefficients  each  of  magnitude  1.  Thus,  for  1p(z)I 
to  be  close  to  constant  on  !zl  =  l,  that  constant  must  be  n1'2.  The 
question  of  the  existence  of  such  polynomials  is  a  classic  one  in 
mathematical  analysis.  Its  study  was  apparently  initiated  by  Hardy 
Cll,p.  199],  and  furthered  by  Littlewood  C8,9],  Erdos  C5],  Newman 
C 1,2,3, 10]  and  others.  A  basic  result  concerning  these  problems  was 
obtained  by  the  author  C4],  which  paved  the  way  for  solutions,  by 
Korner  C7]  and  Kahane  C6],  of  two  of  the  fundamental  conjectures  in 
this  area.  Ue  modify  the  construction  given  in  [4]  to  produce  nearly 
ideals  filters  with  one  notch. 


II.  APPROACH  AND  RESULTS 


Our  starting  point  is  the  polynomial  P,  of  degree  N2-l,  given 


P  P 

L  sx 


p  (  2TTi  j  kN“ 1 )  z-'-*"'*,  :=exp  <2T7iv). 


kzO  j  -0 


It  is  shoun  in  C4D  that  P(^)  satisfies: 


(i)  |P(jN_a)|=N  for  all  integers  j; 

(ii)  For  any  N-1<<<l/2,  i  P  ( )  i  =N+E  for  <<0<l/2,  where 
lEl  <1+277-^5  (  no-1 ; 

(iii)  min  lp(6)i=0(l);  and 

(iv)  1  P<e>  l  <  (2-f37T-‘  )N  +  0(1)  for  all  e. 


Recent  numerical  evidence  suggests  that  (iv)  can  be  strengthened  to: 


Civ’)  i  P O)  i  <  1 . 3N  for  all 


Figures  1  and  2,  which  show  iP(-d)i  as  a  function  of  for  N=10 
(i.e.,  P  of  degree  99)  and  N=60  (P  of  degree  3599),  clarify  the  above 
properties.  Thus,  as  the  magnitude  of  P  is  asymptotically  close  to 
constant  except  for  the  immediate  neighborhood  of  one  point.  By  a 
simple  change  of  variables) it  is  clear  that  this  special  point  can  be 
taken  anywhere  on  the  unit  circle. 


If  P(O)  is  changed  by  removing  the  first  N  terms,  all  of  whose 
coefficients  are  +1,  and  then  dividing  by  there  results 


Q<e> 


N-/  U-i 

■r  e 

j=° 


»xp ( 2TT  i  jkN-1 )  z  j-«  •*-»  >~. 


This  Q,  which  is  of  degree  Na-N-1,  is  the  desired  modification  of  P.  In 
fact,  estimates  (ii)  and  (iv’)  remain  true  for  Q,  while  in  addition 


Q  ( 0 ) =0 


'■j  *V*  «  »  >irx,Hiw  K”  *r  *  -  k~  v  v  «.*  mip.-jsr.’vjn 


Q ( 0 ) —0 .  Also,  it  can  be  shown  that  the  null  width  of  Q  is  less  than 
2/N.  Figures  3  and  4,  which  exhibit  lQ(£)i  as  a  function  of  -*4  for  N=10 
(Q  of  degree  89)  and  N=60  (Q  of  degree  3539>,  and  figure  5,  which 
transforms  the  plot  of  figure  4  to  a  dB  scale,  show  that  Q  is  indeed  the 
nearly  ideal  notch  filter  discussed  earlier.  Once  again,  it  is  clear 
that  a  change  of  variables  allows  the  relocation  of  the  notch  to  any 
desired  0. 
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